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MTI  SYSTEM  SIMULATION  AND 
CLUTTER  OUTPUT  STATISTICS 


INTRODUCTION 

A  radar  return  from  a  patch  of  clutter  usually  consists  of  a  large  number  of  echoes  from  indi¬ 
vidual  scatterers.  Each  of  these  scatterers  moves  randomly  and  introduces  a  randomly  distributed 
doppler  spectrum.  Therefore,  to  simulate  the  performance  of  a  MTI  system,  we  must  sum  a  large 
number  of  randomly  distributed  samples  from  each  radar  pulse  return.  The  computer-time  required  in 
this  case  is  very  lengthy,  and  therefore  it  is  not  feasible  for  us  to  use  this  approach  in  a  large-scale 
simulation.  In  this  report  we  propose  to  use  a  simplified  version  of  the  clutter  model.  Such  a  model, 
as  we  will  show,  will  generate  a  correlation  function  and  a  spectral  density  function  which  are  experi¬ 
mentally  measured.  Such  a  model  reduces  the  required  computer-time  a  great  deal. 

The  second  problem  we  will  investigate  is  the  improvement  factor  of  the  MTI  system.  This 
improvement  factor  is  a  function  of  the  clutter  output  and  will  be  used  as  a  base  for  the  design  of  an 
MTI  system  and  its  required  performance.  The  clutter  output  is  the  weighted  summation  of  a  number 
of  delayed  radar  returns.  Since  the  spectral  distribution  of  these  radar  returns  is  random,  the  improve¬ 
ment  factor  is  also  a  random  function.  In  the  past,  the  expected  value  of  this  improvement  factor  was 
used  for  the  radar  MTI  design;  therefore  it  is  possible  that  it  accounts  for  the  fact  that  a  well-designed 
MTI  system  may  not  be  adequate  to  eliminate  clutter  noises  at  all  times.  In  this  report,  we  present,  as 
an  example,  the  probability  distribution  of  the  improvement  factor  of  a  MTI  system.  This  MTI  system 
uses  binomial  weights,  and  the  results  are  obtained  through  computer  simulation. 

CLUTTER  SPECTRAL  DENSITY  FUNCTION 

The  radar  clutter  returns  usually  are  characterized  by  two  parameters: 

a.  The  radar  cross  section  and 

b.  The  spectral  density  function. 

The  radar  cross  section  is  the  measure  of  the  reflection  of  a  clutter  patch  when  it  is  illuminated  by 
radar  emission.  The  spectral  density  function  is  the  distribution  of  the  random  motion  of  the  scatterers. 
In  general,  since  a  clutter  patch  consists  of  a  large  number  of  scatterers,  and  the  radar  echo  is  the  sum¬ 
mation  of  the  reflections  from  each  individual  scatterer,  the  reflection  from  each  scatterer  is  random 
and  independent.  Therefore,  according  to  the  central  limit  theory,  the  summation  of  these  independent 
random  variables  is  asymptotically  normal  as  the  number  of  random  variable  approaches  infinity. 
Hence,  the  radar  return  changes  into  a  Rayleigh  distribution.  This  fact  has  been  pointed  out  very  early 
by  Lawson  and  Uhlenbeck  til;  later,  other  distributions  such  as  log  normal  have  been  reported  for  sea 
and  land  returns  [2],  The  radar  cross  section  of  clutter  is  a  very  important  factor  in  the  detection  of 
target.  If  the  clutter  is  strong  enough,  it  might  mask  the  target;  however,  some  means  must  be  pro¬ 
vided  to  eliminate  the  clutter  interference.  Since  the  clutter  correlation  time  usually  is  very  long,  the 
integration  scheme  does  not  work  under  this  condition.  To  reduce  clutter  interferences  we  must  avail 
ourselves  of  the  evident  fact  that  the  clutter  in  general  is  stationary.  We  can  use  a  kind  of  velocity 
filter  which  may  filter  out  the  stationary  clutter  and  enhance  the  detection  of  a  moving  target.  This  has 
been  proven  and  successfully  carried  out.  To  design  such  a  velocity  filter,  it  is  very  important  that  we 
possess  knowledge  of  the  characteristics  of  clutter  movement. 
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In  reality,  the  clutter  is  not  completely  stationary.  As  we  pointed  out  earlier,  the  clutter  consists 
of  a  large  number  of  scatterers,  and  since  each  scatterer  moves  randomly  because  of  prevailing  wind  or 
.other  effects,  the  motion  of  the  scatterers  introduces  doppler  shift  of  the  radar  returns.  The  doppler 
shift  is  a  random  function  that  is  characterized  by  the  correlation  function  or  the  spectral  density  func¬ 
tion.  In  this  section,  based  on  the  model  that  a  clutter  patch  consists  of  a  large  number  of  clutter 
*  scatterers  moving  randomly,  we  determine  the  correlation  function. 

Let  us  assume  that  there  is  a  large  number  of  scatterers  located  in  a  patch  of  clutter,  and  let  us 
also  assume  that  point  0  is  the  phase  center  (a  reference  point)  of  this  clutter  patch;  if  a  radar  signal  is 
reflected  from  a  point  scatterer  a,  then  this  reflected  wave  can  be  represented  as  (see  Fig.  1): 

na(t\)  =  ya  cos  («*>o  t  +  2 kR  +  kR  ■  Ti  +  {3a ) 
where  ya  is  the  amplitude  of  the  reflected  wave. 

w0  is  the  radar  carrier  frequency,  for  simplification  we  have  assumed  that  the  radar  signal 
is  narrow  band. 

R  is  the  distance  from  radar  site  to  the  phase  center  point  0. 

R  is  a  unit  vector  of  R. 

r{  is  position  vector  of  the  scatterer  point  a . 

Pa  is  the  phase  due  to  reflection  from  the  scatterer. 
k  is  the  wave  number. 

For  this  information  we  assumed  that  R  »  r\. 


Fig.  1  —  Scattering  by  randomly  moving  clutter  particles 
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At  time  t2  this  particle  moves  to  a  new  position;  therefore,  the  radar  return  becomes 
r\a(t2)  =  ya  cos  (ojq  t2  +  2kR  +  kR  ■  r2  +  (ia). 


Since  r2  =  ri  +  r'  we  have 

Va(h )  =  ya  COS  (w0  t2  +  2 kR  +  kR  ■  (rj  +  r')  +  y3Q). 

Since  kR  represents  a  constant  phase  delay  for  all  scatterers,  it  can  be  ignored.  Expression  -  rt  +  (3a 
is  a  random  phase  which  varies  from  scatterer  to  scatterer;  therefore,  we  can  represent  it  by  a  random 
variable  </>„.  Furthermore,  because  the  time  elapse  from  t\  to  t2  is  very  short,  we  can  assume  that  the 
scatterer  a  moves  at  a  constant  velocity.  Then  kR  •  r'  =  (2i t/\)  v  (t2  -  tx),  and  v  is  the  radial  com¬ 
ponent  of  the  velocity  of  scatterer  a.  We  should  notice  that  v  is  a  random  variable. 


The  correlation  function  is: 


=  E 


ya  y? 

a  $ 


~2  C0S  (W0  ■”  h )  —  (27t/Xv)(^2  “  t\)  +  (<t>a  “  0/?)) 


+  j  COS  (co0  (tl  +  t2)  +  (<t>a  -  <f>l 3)  +  (2tt/Xv)(/2  ~  ^i)) 


The  reflection  coefficient  y,  phase  angle  </>,  and  velocity  v  are  all  independent  random  variables,  and  the 
process  is  stationary.  The  phase  angle  </>  has  a  uniform  distribution  between  ~n  to  7 r,  because  we 
assume  that  the  scatterers  are  randomly  distributed  within  the  clutter  patch.  Because  of  this  assump¬ 
tion,  the  correlation  function  becomes 


where 


R  (r)  —  R0  cos(w0r  —  (27r/X)vr)j 


t  -  h  -  h 

*0=  }  ]T  y\. 

a 

The  velocity  v  may  contain  two  components.  A  constant  component  which  is  common  to  each 
scatterer,  and  which  represents  the  velocity  that  the  clutter  patch  moves  as  a  single  body.  The  second 
component,  which  varies  randomly  from  scatterer  to  scatterer,  represents  the  velocity  of  movement  of 
each  scatterer,  thus 


The  correlation  function  is 


v  =  v0  +  Av. 


R(t)  = 


Where  fQ  is  the  doppler  due  to  the  mean  velocity  v0,  and  A/  is  the  random  component.  Let  us 
suppose  that  this  signal  has  been  converted  into  the  base  band,  and  let  us  assume  that  the  probability 
density  function  A/ has  a  probability  function  S(f),  then 

R  (r)  =  R0  f  S (A/)  cos(2tt(/0  + 
or 

R  (r)  =  Rejtfoe'2^07  f  S  j 
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which  states  that  the  correlation  function  is  the  Fourier  transformation  of  the  S(f).  The  S(f)  function 
is  commonly  referred  to  as  the  power  spectral  density  function  of  the  random  process  r?  (f).  We  may 
also  notice  that  function  S(f)  represents  the  probability  density  function  of  the  doppler  frequency  of 
the  clutter  scatterer.  The  clutter  echo  can  thus  be  simulated  by  the  summation  of  many  scatterer 
returns. 

SIMULATION  OF  CLUTTER  RETURN  SIGNAL 

As  we  have  seen  in  the  previous  section,  we  can  simulate  the  clutter  return  by  summing  a  large 
number  of  scatterer  returns,  each  of  which  has  a  random  doppler  frequency  £(/),  as  its  probability 
density  function.  However,  in  the  actual  simulation,  this  approach  requires  too  much  computer-time 
and  therefore  is  not  feasible  for  a  large-scale  simulation.  In  a  real  radar  system,  most  MTI  or  doppler 
filter  is  performed  digitally.  The  returns  of  the  radar  echo  at  each  range  bin  are  sampled  and  processed. 
Therefore,  if  these  samples  were  used  for  simulation,  the  results  would  represent  an  actual  radar  opera¬ 
tion.  It  is  evident  that  such  an  approach  requires  much  less  computation  because  one  sample  is  all  that 
is  required  instead  of  the  summation  of  a  large  number  of  returns.  Furthermore,  the  measured  clutter 
data  is  usually  characterized  by  the  correlation  function  which  was  computed  from  these  sampled  data  at 
each  range  bin.  A  numerical  Fourier  transformation  of  this  correlation  function  is  then  performed  to 
obtain  the  spectral  density  function  [3].  It  is  evident  that  if  we  can  reconstruct  these  sampled  data  from 
their  computer  spectral  density  function,  these  data  can  be  used  to  accurately  simulate  the  radar  MTI 
operations.  In  this  section  we  shall  show  how  these  sampled  data  can  be  reconstructed  from  the  known 
spectral  density  function. 

The  clutter  return  is  known  to  be  a  wide  sense  stationary  process.  Its  mean  value  is  independent 
of  time,  and  its  correlation  function  is  only  a  function  of  the  correlation  time. 

The  appendix  of  this  report  shows  a  time  function  of  the  following  form  that  can  be  used  to 
represent  such  a  random  process: 

x(f)  =  cos(27r  ft  +  £). 


Both  /  and  £  are  random  variables.  Random  variable  £  has  uniform  distribution  from  —tt  to  7 r, 
while  /  has  a  probability  density  function 

p(f)  -  2  S(f), 

where  S(f)  is  the  spectral  density  function  which  has  a  correlation  function 

R  (r)  =  E[x(t)x(t  +  t)]. 


This  time  function  can,  then,  be  used  as  the  radar  return  for  simulation  purposes.  Frequency  / 
represents  the  doppler  shift,  and  £  represents  the  initial  phase.  We  notice  that  each  x(t)  function  for  a 
given  /  and  £  represents  a  sample  of  the  ensemble  of  the  random  stochastic  function  of  the  clutter 
return.  For  a  different  sample,  the  /  and  £  vary  according  to  each  probability  density  function,  as 
shown  above. 


MTI  IMPROVEMENT  FACTOR 


An  MTI  system,  as  shown  in  Fig.  2,  has  a  filter  weight  an  and  a  delay  of  r  at  each  stage.  The 
clutter  output  of  such  a  system  is  equal 


So  = 


£  a„  exp  (jlirfiM )  2 
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Fig.  2  —  An  MTI  system 


Where  ft  is  the  doppler  frequency  of  the  target,  r  is  the  delay,  and  it  is  assumed  that  the  input  signal 
has  a  unity  amplitude.  We  notice  that  the  pulse  repetition  frequency  (PRF)  fr  is  equal  to  1/t.  If  the 
target  doppler  is  an  integer  of  /r,  then  the  input  signal  function  S0  repeats  itself.  Sq  is  thus  a  periodic 
function.  For  convenience,  we  may  introduce  a  normalized  doppler  frequency  fh  that 

fl-fjfn 


then  the  signal  output  becomes 

o  _  ^  ^  r, 

n  m 

When  the  target  velocity  is  not  known  a  priori,  we  may  assume  that  the  target  velocity  has  a  uniform 
distribution.  Under  this  condition,  the  average  signal  gain  is  thus 

£  (S0)  =  5>„2. 

n 

Similarly  the  clutter  output  is 

/o  tt\)  f t 

C  =  anam  e 

n  m 

where  C0  is  amplitude  of  the  input  clutter.  If  the  clutter  has  a  spectral  density  function  Gif),  then  the 
average  clutter  output  is 

EiC)  =  C02  ££  anamRnm, 

n  m 

where  Rnm  is  the  correlation  function  and  it  is  the  Fourier  transformation  of  the  spectral  density  func¬ 
tion. 


The  improvement  factor  is  defined  as 


LF=  E 


Input  target  signal  power/input  clutter  power 
Output  target  signal  power/output  clutter  power 


n 
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If  the  filter  weight  is  normalized  so  that 


2>„2=i 

n 

then 


nm  • 


The  correlation  function,  Rnm,  is  the  Fourier  transformation  of  the  spectral  density  function. 
According  to  the  data  measured,  most  clutter  spectral  density  has  a  normal  distribution.  Although 
some  ground  clutter  has  been  reported  to  have  a  f  distribution,  it  has  been  shown  that  it  has  a  general 
bell  shape  as  that  of  a  normal  distribution,  except  for  a  long  tail.  We  shall  assume,  however,  that 
clutter  spectral  has  a  Gaussian  distribution,  and  furthermore,  that  the  clutter  mean  velocity,  because  of 
prevailing  winds  or  platform  movement,  usually  has  to  be  removed  by  some  other  means.  In  our  fol¬ 
lowing  discussion,  this  mean  velocity  is  assumed  to  be  zero.  The  clutter  spectrum  density  function  has 
the  following  form 


G(f)  = 


V2 


7T 


-r 

la2 


For  convenience,  the  variable  /  and  standard  deviation  cr  are  both  normalized  with  respect  to 
PRF.  We  therefore,  can  consider  a  range  of /from  0  to  1.  The  correlation  function  has  the  form 

Rnm  =  exp[—2n2cr2(n  —  m )2]. 

A  set  of  improvement  factors  for  MTI  filter  with  binomial  weights  is  shown  in  Fig.  4.  The  MTI 
improvement  is  plotted  in  dB  scale  and  it  is  plotted  as  a  function  of  the  normalized  cr  for  a  value  less 
than  0.25,  or  a  standard  deviation  of  doppler  no  more  than  1/4  of  the  PRF. 

The  binomial  coefficients  are  used  for  MTI  filter  weights.  It  has  been  shown  that  binomial  weight 
is  only  optimal  at  cr  =  0  [4].  Nevertheless,  it  has  been  used  since  the  beginning  of  the  MTI  system, 
and  it  is  very  convenient  to  use.  There  are  other  types  of  filter  weights  [4-6]. 

We  have  noticed  in  the  previous  formulation,  that  the  target  velocity  is  assumed  to  be  uniformly 
distributed;  therefore,  if  the  target  velocity  is  known,  we  might  use  some  other  filter  weights  to  achieve 
a  better  signal  gain.  One  such  example  is  the  doppler  filter  [6], 

MTI  SIMULATION 

We  assume  that  the  clutter  return  signal  for  the  MTI  simulation  has  the  following  form: 

/  —  channel  x(t)  —  cos(27t ft  +  f ) 

Q  —  channel  x(t)  =  sin(27r/f  4- 

where  /  and  £  are  random  variables.  Random  variable  /  has  a  normal  distribution  with  zero  mean  and 
a  given  cr.  While  £  has  a  uniform  distribution  from  —it  to  tt.  The  MTI  output  is  then 

I  —  channel 

Q  =  £  cos (nfr  +  |) 


Q  —  channel 

CQ  =  £  an  sin  (nfr  +  0 
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where  r  is  interpulse  delay  time.  The  total  clutter  output  is  then 
C2  =  C,2  +  Cl 

=  ££  anam  cos (nft  +  £)  cosimft  +  o 

n  m 

+  ££  anam  sin(w/,f  +  0  smimft  +  £) 

n  m 

=  ~  ££  a„am{[cos((n  -  m)/r)  -  cos ((m  +  n)/r  +  2^)] 

n  m 

+  [cos  ((«— /?z)/t)  +  cos  ((m  +  «)/t  +  2^)]} 

=  anam  cos ((«  -  m)/r). 

n  m 

If  the  weights  an' s  are  normalized,  then  the  expected  value  of  the  reciprocal  represents  the  improve¬ 
ment  factor. 

Figure  3  shows  the  results  of  one  of  such  simulation.  This  figure  shows  that  we  have  plotted  the 
improvement  factor  as  a  function  of  the  normalized  cr,  the  variance  of  the  spectrum  density  function. 
Also,  this  figure  shows  that  we  have  plotted  the  expected  value  of  the  improvement  factor  of  a  three- 
pulse  MTI  system.  Later  at  each  <x  value,  we  have  computed  a  large  sample  of  clutter  output.  Each 
sample  has  a  different  doppler  frequency  which  is  generated  from  a  normal  distribution  with  a  variance 
of  cr.  The  initial  phase  £  was  selected  from  a  uniform  distribution  —7 r  to  tt .  We  computed  the 
improvement  factor  of  a  total  of  1000  samples  for  each  cr  in  Fig.  3,  and  then  we  computed  its  average 
value  which  is  also  plotted  in  Fig.  3.  We  may  see  that  the  average  improvement  factor  value  is  almost 
identical  to  that  of  the  theoretical  value.  This  may  prove  the  validity  of  this  simulation  approach.  Fig¬ 
ure  4  shows  the  average  improvement  factors  as  a  function  of  normalized  cr  for  2  to  6  pulse  MTI  sys¬ 
tems. 


BD. 


Fig.  3  —  Improvement  factor  vs  normalized  clutter  a- 
for  theoretical  and  simulation  results 
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Fig.  4  —  Average  improvement  factor  as  a  function  of 
clutter  spectral  standard  deviation 


STATISTICAL  DISTRIBUTION  OF  THE  IMPROVEMENT  FACTOR 

In  the  previous  section  we  demonstrated  how  the  improvement  factor  is  computed.  We  also 
showed  how  the  simulation  of  improvement  factor  is  formed.  This  improvement  factor  usually  is  used 
as  a  base  for  the  design  of  an  MTI  system.  However,  this  improvement  factor  is  a  random  function. 
In  the  past  its  expected  value  has  been  used  for  this  purpose;  however,  in  a  real  system,  this  value  may 
fluctuate.  There  are  two  reasons  that  contribute  to  the  fluctuation.  First,  in  computing  the  improve¬ 
ment  factor  we  used  the  average  target  velocity.  In  reality  however,  the  target  velocity  may  fluctuate 
therefore  producing  a  fluctuation  in  achievable  improvement  factors.  Second,  the  clutter  doppler  is  a 
random  variable;  hence  the  clutter  output  is  also  a  random  function.  It  varies  from  clutter  sample  to 
clutter  sample  and  it  also  varies  at  different  sampling  time  (if  sampling  interval  is  greater  than  correla¬ 
tion  time). 

The  variation  of  target  velocity  is  not  addressed  here.  In  general,  this  problem  can  be  easily 
solved  by  using  some  sort  of  doppler  filter. 

In  this  section  we  address  the  second  problem.  The  clutter  output  is  a  summation  of  random 
variables.  Usually  we  can  assume  that  such  a  function  has  a  normal  distribution  if  the  summed  random 
variables  are  independent.  Unfortunately,  the  clutter  output  is  the  summation  of  many  MTI  pulse 
returns  which,  because  of  long  clutter  correlation  time,  cannot  be  treated  as  independent.  The  assump¬ 
tion  of  a  normal  distribution  is  therefore  invalid.  There  is  no  known  probability  density  function  to 
describe  such  a  process;  therefore,  the  following  results  are  based  on  the  computer  simulations. 

We  have  performed  this  simulation  for  a  3-pulse  MTI  system  with  binomial  weights.  We  vary  the 
normalized  cr  from  0.01  to  0.1.  For  each  cr,  we  compute  10,000  samples.  For  each  sample,  we  gen¬ 
erate  the  doppler  frequency  /  and  phase  f .  We  generate  the  doppler  frequency  according  to  a  normal 
density  function  with  zero  mean  and  a2  as  the  variance.  Then  we  generate  the  random  number  £  as  a 
uniform  distribution  from  ~i r  to  it.  Finally  we  compute  the  clutter  output.  This  computed  clutter  out¬ 
put  is  then  normalized  with  respect  to  the  average  clutter  output.  We  plot  the  cumulative  probability  of 
such  samples  as  shown  in  Fig.  5.  There  are  five  such  curves,  each  for  a  given  c r  value  which  varies 
from  0.01  to  0.1.  One  of  the  interesting  results  that  we  may  notice  is  that  these  cumulative  probability 
curves  are  not  a  function  of  the  normalized  cr.  Similar  results  were  obtained  for  4-pulse  and  5-pulse 
MTI  system.  The  clutter  level  in  this  plot  is  normalized  with  respect  to  the  average  clutter  level.  In  a 
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Fig.  5  —  Cumulative  probability  of  normalized  clutter  output  level  (nor¬ 
malized  with  respect  to  average  clutter  output)  for  a  3-pulse  MTI  system 


three-pulse  canceller,  as  shown  in  Fig.  4,  the  average  clutter  output  amplitude  at  cr  =  0.05  is  0.067  (or 
—23.48  dB).  In  Fig.  5  we  see  that  80%  of  the  samples  have  less  than  this  clutter  output.  However,  if 
we  wish  to  keep  a  90%  probability  improvement  factor,  the  allowable  clutter  becomes  0.134  (or  —17.46 
dB). 


Figure  6  shows  the  same  cumulative  probability  curve  for  2-,  3-,  4-,  5-,  and  6-pulse  MTI  systems. 
This  figure  shows  the  same  normalized  clutter  as  that  shown  in  Fig.  5.  For  example,  for  n  =  2,  the 
probability  that  the  clutter  output  is  less  than  or  equal  to  the  average  clutter  is  0.66.  Table  1  lists  the 
probabilities  of  MTI  filter  which  are  less  than  or  equal  to  its  expected  values  for  2  to  6  pulses. 
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Fig.  6  —  Cumulative  probability  of  normalized  clutter  output  level 
(normalized  with  respect  to  average  clutter  output) 
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Table  1  —  Probabilities  of 
MTI  Filter 


No.  of  Pulses 

Probability 

2 

0.66 

3 

0.80 

4 

0.87 

5 

0.92 

6 

0.95 

It  is  evident  that  the  increase  in  the  number  of  pulses  not  only  reduces  the  average  clutter  output,  but 
there  is  also  a  higher  probability  that  the  clutter  output  will  not  exceed  that  level. 

The  information  contained  in  Fig.  6,  and  in  Fig.  4,  can  be  used  to  find  the  probabilistic  improve¬ 
ment  factor  for  a  conventional  binomial  weighted  MTI  system.  For  example,  at  a  normalized  a  =  0.05, 
the  average  improvement  factor  of  a  4-pulse  canceller  is  32  dB.  From  the  information  in  Fig.  6  we  may 
notice  that  0.87  probability  of  all  samples  will  achieve  this  improvement  factor.  However,  if  we  want  to 
achieve  a  0.95  probability,  the  improvement  factor  decreases  to  26.4  dB. 

CONCLUSION 

In  this  report  we  have  reviewed  how  a  clutter  correlation  function  is  formulated  and  how  it  relates 
to  the  spectral  density  function.  To  simplify  MTI  simulation,  we  have  formulated  a  random  time  func¬ 
tion  that  can  produce  a  correlation  function  and  spectral  density  function  as  the  one  measured  experi¬ 
mentally.  We  have  checked  the  validity  of  such  simulation  against  the  theoretical  results;  we  have 
computed  the  improvement  factors  for  binomial  weighted  MTI  systems  and  have  plotted  its  cumulative 
probabilities  of  staying  below  a  certain  clutter  output  level. 
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APPENDIX 


Given  a  random  function  x(r),  its  expected  function  is 
£[*(/)]  =  £ [cos (27r ft  T-  £ )  ] 

=  £[005277/1?  cos  £  —  sin  27r  ft  sin  £]. 

Since  /  and  f  are  independent, 

E[x{t)]  =  £(cos277/f)£(cos4:)  —  Eisinlnft  E( sin£). 

If  x{t)  is  stationary,  £[x(r)]  must  be  a  constant  and  not  a  functin  of  t.  This  will  happen  when 

E( cos£)  -  £(sin£)  =  0. 

This  implies  that  the  characteristic  function  <f>(k)  of  £ 

0(1)  =  0, 

where 

<f>ik)  =  f  p(0  eJk*d( 
and  p{£)  is  the  probability  density  function  of  random  variable 


The  correlation  of  x  it)  is 

R  (r)  -  £  (cos  (277 /r  +  £)  cos(2t7  fit  +  r)  +  £) 

=  y  £(c0s(277/t))  +  ~  E  (cOs(477  ft  +  277 fr  +  2£)) 

if  £(cos2£)  =  £(sin2£)  =  0. 


Then 


R  (r)  —  y  £(cos2t7/t) 
or 

£  (t)  =  y  J*  j?(/)  cos277  frdf, 

where  p(f)  is  the  probability  density  function  of  the  frequency  /.  Since  R  (r)  is  the  Fourier  transfor 
mation  of  the  spectrum  density  function  Sit),  we  have 

Si/)=jp(f). 


The  probability  density  function  of  £  must  have  a  uniform  distribution  in  the  range  from  —77  to  77 
This  can  be  shown  as  follows: 


4>ik) 


<f>(k) 


JL  J*, 

277 


sin/:  77 

kn 


0.  k  =  1,  2. 
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